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Abstract-  A  finite  element  muscle  model  to  study  the  principles  of 
intermuscular  myofascial  force  transmission  is  developed.  The 
results  obtained  explain  force  differences  at  the  distal  and 
proximal  tendons  of  muscles  that  have  mechanical  interaction. 
This  is  in  agreement  with  experimental  findings  in  other  recent 
studies.  The  strain  distributions  found  along  the  fiber  direction 
indicate  intermuscular  myofascial  force  transmission.  A 
consequence  is  that  active  force  generated  within  one  muscle  may 
be  exerted  at  the  tendon  of  another  muscle. 
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I.  INTRODUCTION 


Performing  the  complex  task  of  controlled  bodily 
movements  requires  the  transmission  of  the  force  generated  by 
the  skeletal  muscles  onto  the  skeleton.  Although  for  this 
transmission,  the  myotendinous  junction  is  widely  accepted  to 
be  the  major  site,  rather  recently  the  transmission  from  the 
muscle  fibers  onto  the  intramuscular  connective  tissue  has  been 
proven  to  be  important  as  well  [1-5].  Referring  to  this  kind  of 
transmission  as  myofascial  force  transmission,  Huijing  [6] 
suggested  that  force  could  be  transmitted  from  the  muscle  by 
this  path  as  well  (extramuscular  myofascial  force 
transmission).  Recent  experiments  revealed  significant 
changes  in  the  length-force  characteristics  of  rat  extensor 
digitorium  longus  (EDL)  muscle  on  interference  with  [7]  or 
interactions  with  the  surrounding  extramuscular  tissues.  It  was 
also  shown  that  adjacent  muscles  interact  mechanically  [8].  In 
the  present  work,  the  two-domain  finite  element  muscle  model 
described  recently  [9]  is  extended  to  include  the  principles  of 
intermuscular  myofascial  force  transmission.  In  experimental 
work  on  the  subject,  the  muscle  is  not  accessible  for  local 
strain  measurements,  as  it  has  to  remain  surrounded  by 
compartmental  connective  tissues.  This  modeling  study  aims  at 
providing  such  analysis  of  local  strain  to  enhance 
understanding  of  this  concept. 

II.  METHODOLOGY 


A  3D-finite  element  muscle  model  (linked  fiber-matrix 
mesh  model:  lfmrn  model)  with  a  two  domain  approach  was 
developed  [9].  This  model  consists  of  two  meshes  that  are 
linked  elastically  representing  the  extracellular  matrix  (matrix 
mesh)  and  muscle  fiber  domains  (fiber  mesh),  which  occupy 
the  same  space.  These  two  meshes  are  built  using  earlier 
developed  myofiber  or  extracellular  matrix  elements  that  are 
introduced  into  the  finite  element  program  ANSYS  5.5.1  as 
user  defined  elements.  The  two  meshes  are  rigidly  connected 
to  single  layers  of  elements  representing  the  muscles’ 
aponeurosis,  for  which  a  standard  element,  HYPER58  from  the 

element  library  of  ANSYS  5.5.1  is  used. 
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Fig.l.  The  deformed  shape  of  the  model.  The  model  consists  of  two  EDL 
models.  One  that  is  restrained  at  a  low  length  (i.e.  2.5  mm  shorter  than  the 
initial  length),  which  has  interaction  with  another  EDL  model  that  is 
lengthened  distally  (from  Aim  =  -2.5mm  to  Aim  =  1.5  mm).  Muscular 
geometry  at  the  initial  muscle  length  (lm  =  28.7  mm)  is  represented  by  dotted 
lines.  A  3D  local  coordinate  system  is  used  for  the  analysis  and  presentation  of 
the  results. 


In  the  present  study,  using  the  same  methods,  a  model 
representing  whole  EDL  isolated  from  the  surrounding 
connective  tissues  was  built.  As  a  subsequent  step,  two  of  such 
models  were  connected:  the  corresponding  nodes  of  the  matrix 
meshes  of  the  two  models  were  linked  elastically.  Such  links 
represent  the  intermuscular  connections,  providing  a  means  to 
show  the  principles  of  mechanical  interaction  between  the  two 
muscles  and  a  pathway  for  intermuscular  force  transmission. 

To  assess  the  principles  of  such  interaction  and  force 
transmission,  one  EDL  model  referred  to  as  restrained  muscle 
was  kept  at  a  low  length  ( Aim  =  2.5  mm  below  the  original 
length),  while  the  other  model  (referred  to  as  lengthened 
muscle)  was  lengthened  from  this  low  length  up  to  1.5  mm 
over  the  original  length  (Fig.  1).  Throughout  the  analysis,  both 
modeled  muscles  were  maintained  maximally  activated.  The 
length-force  characteristics  of  the  isolated  model  will  be 
compared  to  experimental  data  for  EDL  in  such  conditions. 
Stress  and  strain  distributions  of  the  isolated  model  and  the 
model  with  intermuscular  links  will  be  compared  as  well. 

III.  RESULTS 

Isolated  EDL  Model 

The  length-  total  force  characteristics  of  modeled  isolated 
EDL  and  experimental  data  [7]  are  shown  in  Fig.  2.  The  model 
forces  are  maximally  5%  above  or  below  the  experimental 
ones  suggesting  a  fairly  good  agreement.  For  the  highest 
lengths  considered,  the  model  forces  remain  high  whereas  the 
experimental  force  decreases  with  increasing  length.  This 
could  be  due  to  over  estimated  passive  length-force 
characteristics  of  the  model  as  well  as  the,  unrealistic  modeling 
assumption  of  an  initially  uniform  sarcomere  distribution. 
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Fig.2.  The  isometric  muscle  length-total  force  curve  for  isolated  rat  EDL 
muscle.  Experimental  EDL  data  is  adopted  from  [7].  The  muscle  force  (Fm) 
is  normalized  for  optimum  force  (Fmo).  Initial  muscle  length  equals  28.7  mm. 

Model  with  Intermuscular  Connections 

The  muscle  length-  total  force  characteristics  of  the 
lengthened  and  the  restrained  muscle,  which  have 
intermuscular  links  between  their  matrix  meshes  representing 
intramuscular  collagen  is  presented  in  Fig.  3.  Note  that  for 
distributions  of  strain  and  stress  in  the  fiber  mesh  of  the 
restrained  muscle  are  provided  for  three  lengths  of  the 
lengthened  muscle  (Fig.  4).  At  Aim  =  -2.5  mm  (  i.e.  both 


Muscle  Force  (Fm/Fmo) 


Muscle  Length  (mm) 


b 
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Fig.3.  Comparison  of  proximal  and  distal  forces  of  the  modeled  EDL  muscles. 

(a)  Length-total  force  curves  for  the  lengthened  muscle.  The  muscle  force 
(Fm)  is  normalized  with  respect  to  optimum  distal  force  (Fmo). 

(b)  Total  force  for  the  restrained  muscle  as  a  function  of  the  lengthened 
muscles’  length.  The  muscle  force  is  normalized  for  optimum  proximal  force. 
Initial  muscle  length  equals  28.7  mm 


muscles  are  at  the  lowest  length),  the  fiber  direction  strain 
( E-n  )  does  not  exhibit  variations,  i.e.  the  lengths  of  sarcomeres 
in  series  within  fibers  are  highly  uniform  (Fig.  4a).  As  the 
lengthened  muscle  reaches  higher  lengths  ( Aim  =  +1.5  mm) 
this  homogeneity  is  replaced  by  on  average  10%  longer 
sarcomeres  at  the  distal  end  of  muscle  fibers  than  at  their 
proximal  end. 

The  consequences  of  these  differences  in  strain  are  clearly 
visible  on  the  fiber  direction  stresses  (  CU2  )•  If  both  muscles 
are  at  low  length,  the  stresses  along  the  fibers  of  the  restrained 
muscle  are  rather  uniform.  As  the  adjacent  muscle  is 
lengthened,  distributions  of  (J,2  become  apparent:  at  distal 
ends  of  fibers  approximately  50%  higher  stresses  are 
encounterd  than  their  proximal  ends. 

For  a  further  assessment  of  mechanics  of  interaction 
between  the  muscles,  the  fiber  direction  strains  and  stresses 
within  the  fiber  mesh  at  the  highest  length  studied  (i.e.  bn  = 
+1.5  mm)  are  considered  for  the  lengthened  muscle  of  the 
modeled  pair,  in  comparisson  to  the  isolated  EDL  model  (Fig. 
5).  A  clear  difference  in  strain  distributions  between  the  two 
muscles  is  visible  (Fig.  5a).  Except  for  the  minor  distributions 
in  the  most  proximal  part  of  the  isolated  muscle  (which  is 
ascribable  to  the  asymmetric  EDL  muscle  geometry),  the 
strains  are  highly  uniform  over  the  whole  fiber  mesh.  In 
contrast,  for  the  lengthened  muscle  with  intermuscular 
connections  to  the  restrained  muscle,  a  high  level  of  variation 
in  e21  distribution  is  seen:  The  sarcomeres  along  the  proximal 
ends  of  the  muscle  fibres  are  approximately  20% 


Restrained  Muscle 
a 


b 


Fig. 4.  Distributions  of  fiber  direction  strain  and  stress  for  the  fiber  mesh  of  the 
restrained  muscle  as  a  function  of  the  length  of  the  lengthened  muscle. 

(a)  Distributions  of  fiber  strain  at  three  muscle  lengths. 

(b)  Distributions  of  fiber  stress  at  three  muscle  lengths. 

The  dotted  line  contour  indicates  muscle  geometry  at  the  initial  length. 
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Fig.5.  Comparison  of  strain  and  stress  distributions  for  the  fiber  mesh  of  the 
lengthened  and  isolated  muscles. 

(a)  Strain  distributions  of  the  fiber  mesh  for  the  fiber  (22)  direction. 

(b)  Stress  distributions  of  the  fiber  mesh  for  the  fiber  (22)  direction. 

Both  muscles  are  represented  at  the  high  length  ( Aim  =  +1.5mm).  The  dotted 
line  contour  indicates  muscle  geometry  at  the  initial  length. 


shortened,  while  towards  the  distal  ends,  their  lengths 
gradually  increase,  such  that  the  sarcomeres  at  most  distal  end 
are  up  to  67%  longer.  The  differences  between  the  muscles  are 
also  remarkable  for  the  fiber  direction  stresses  (Fig.  5b). In 
agreement  with  the  strain  distributions  for  the  isolated  muscle, 

the  (5 22  values  are  highly  uniform  (i.e.  approximating  unity  all 
over  the  fiber  mesh).  However,  the  lengthened  muscle  of  the 
pair  exhibits  a  sizable  variation  in  022  distribution,  ranging 
from  aproximately  0.80  to  0.50.  In  that  muscle  stress 
decreases  from  the  proximal  end  of  the  muscle  fibres  towards 
their  distal  ends.  Also  note  that  the  mean  stress  for  the 
isolated  muscle  (which  is  approximately  unity)  is  significantly 
higher  than  that  of  the  lengthened  muscle  with  intermuscular 
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Fig.6.  Comparison  of  length-total  distal  force  curves  of  modeled  muscle  with 
intermuscular  connections  and  isolated  muscles. 

Muscle  force  (Fm)  is  normalized  for  optimum  force  (Fmo)  of  the  muscle  with 
intermuscular  connections. 


connections.  Although  <3 ,  ,  does  attain  values  arround  unity 
at  certain  limited  locations,  stress  values  over  the  major  part  of 
the  fiber  mesh  are  much  lower.  As  fiber  direction  stress  in  the 
fiber  mesh  is  the  dominating  factor  determining  muscle  force, 
this  result  is  highly  interesting  since  the  distal  forces  of  the 
lengthened  muscle  are  higher  than  that  of  the  isolated  muscle 
throughout  the  length  range  and  at  this  particular  length  (Fig. 
6).  Such  higher  forces  in  spite  of  the  lower  mean  stresses  can 
only  be  explained  by  intermuscular  myofascial  force 
transmission  from  the  restrained  muscle  to  the  lengthened  one. 

IV.  DISCUSSION 

Connective  tissue  network  and  force  transmission 

Muscle  can  be  represented  as  an  extensive  3D  set  of 
organized  endomysial  tunnels,  within  which  the  myofibers 
operate  [3].  The  identification  of  intramuscular  force 
transmission  from  the  lateral  perimeter  surface  of  myofibers 
onto  the  connective  tissue  network  (myofascial  force 
transmission)  is  relatively  recent  and  limited.  However, 
various  studies  indicated  [1,3]  or  experimentally  confirmed  [1, 
3-6,  10]  such  force  transmission  and  suggested  an  important 
functional  significance  of  the  connective  tissue  network.  This 
endomysial  network  is  continuous  with  the  perimysium  system 
and  can  be  regarded  as  a  continuous  system  of  tunnels  or  tubes 
as  well.  In  addition,  the  muscular  epimysium  is  continuous 
with  this  system  as  well  as  with  the  connective  tissue  of 
adjacent  muscles  in  addition  to  other  elements  of  the 
compartment  in  which  the  muscles  operate  [e.g.  5,  7]. 
Therefore,  in  vivo  there  is  an  additional  route  for  force 
transmission  out  off  the  muscle,  which  completely  bypasses 
the  tendon  of  the  muscle  that  was  the  source  of  the  generated 
force  [5,  7,  11].  In  the  majority  of  the  earlier  work  on  muscle 
functioning,  muscles  are  investigated  isolated  from  their 
surrounding  [e.g.  1,  12,  13,  14].  Such  a  type  of  experiment 
eliminates  any  possible  mechanism  of  inter  or  extra-muscular 
force  transmission  from  or  to  the  muscle  studied. 

Results  of  the  model 

The  model  developed  in  this  study  provides  the  possibility 
to  interpret  effects  of  myofascial  force  transmission  between 
two  maximally  activated  muscles  with  a  set  of  elastic  links 
between  their  matrix  meshes  representing  their  intramuscular 
connections.  A  finite  element  muscle  model  with  such  an 
approach  is  novel.  The  remarkable  result  of  the  analysis  is  the 
significant  differences  between  the  distal  and  proximal  forces 
of  both  of  the  restrained  and  lengthened  muscles.  These 
differences  are  in  such  that  for  almost  all  of  the  length  range 
considered,  the  distal  forces  of  the  lengthened  muscle  are 
higher  than  the  proximal  forces.  This  proximo-distal  force 
difference  increases  with  increasing  muscle  length.  For  the 
restrained  muscle,  the  distal  force  decreases  as  the  lengthened 
muscle  reaches  higher  lengths.  These  results  are  highly 
compatible  with  the  very  recent  findings  in  experimental  work 
exploring  intermuscular  interaction  via  myofascial  force 
transmission  within  rat  hindlimb.  In  their  study  investigating 
the  effects  of  blunt  dissection  and  compartmental  fasciatomy 
on  rat  EDL  length-force  characteristics,  Huijing  and  Baan  [7] 
showed  that  a  systematic  manipulation  of  the  extra  and  inter¬ 
muscular  connective  tissue  structures  altered  proximal  EDL 
force  remarkably.  These  changes  were  decreases  in  force  and 


4  of  4 


Proceedings  -  23rd  Annual  Conference  -  IEEE/EMBS  Oct.25-28,  2001,  Istanbul,  TURKEY 


changes  in  the  muscle  optimum  length,  as  a  result  of  decreased 
inter  and  extra-muscular  interactions  as  EDL  was  dissected 
further.  This  indicates  the  importance  of  extramuscular 
connective  tissue  in  force  transmission. 

The  advantage  of  the  present  model  is  its  capability  to 
examine  the  mechanics  by  analysis  of  local  stress  and  strain 
distributions.  Such  distributions  reveal  that  even  though  the 
restrained  muscle  is  fixed  at  a  constant  low  length,  the  fiber 
mesh  of  the  restrained  muscle  exhibits  strain  variations  along 
the  fiber  direction  as  the  second  muscle  is  lengthened  (Fig.  4). 
This  is  reflected  on  the  fiber  direction  stresses  of  the  same 
muscle. 

The  model  results  show  that  such  distributions  are  also 
present  in  the  lengthened  muscle  of  the  adjacent  pair  (Fig.  5), 
while  the  isolated  muscle  model  exhibits  highly  uniform  stress 
and  strain  distributions.  As  mean  fiber  direction  stress  within 
the  fiber  mesh  is  considerably  lower  for  the  lengthened  muscle 
than  for  the  isolated  muscle,  the  higher  distal  forces  associated 
with  the  lengthened  muscle  indicate  that  active  force  generated 
within  one  muscle  may  be  exerted  at  the  tendon  of  another 
muscle  by  intermuscular  myofascial  force  transmission.  This 
tendency  of  ‘focusing’  of  active  forces  of  adjacent  muscles  on 
the  distal  tendon  of  a  distally  lengthened  muscle  is  a 
remarkable  conclusion  for  which  further  experimental 
confirmation  will  be  sought  in  the  anterior  tibial  compartment 
muscles  of  rat  hind  limb. 

V.  CONCLUSION 

Present  application  of  the  linked  fiber  matrix  mesh  model 
reveals  the  principles  of  the  significant  role  of  intermuscular 
myofascial  force  transmission  on  the  mechanical  behavior  of  a 
muscle  in  an  in  vivo  situation.  In  agreement  with  experimental 
findings,  it  clearly  demonstrates  the  altered  length  force 
characteristics  with  respect  to  the  isolated  situation  and 
increases  the  understanding  of  the  in  vivo  mechanics.  The 
conclusion  that  active  force  generated  within  one  muscle  may 
be  exerted  at  the  tendon  of  another  muscle  by  intermuscular 
myofascial  force  transmission  is  one  of  the  remarkable  results 
of  this  analysis. 
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